Study of the helix structure of QCD string 
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The properties of a helix-like shaped QCD string are studied in the context of the Lund fragmentation model, 
where the concept of a 3-dimensional structure of the gluon field offers an alternative aproach to the modelling 
of the transverse momentum of hadrons. The paper is focused on the phenomenology of the model, which in- 
troduces correlations between transverse and longitudinal components of hadrons, as well as azimuthal ordering 
of hadrons along the string. The similarities between 2-particle correlations stemming from the helix-string 
structure and those commonly attributed to the Bose-Einstein interference are pointed out. It is shown that the 
charged assymetry observed in hadron production of close hadron pairs can be associated with fluctuations in 
the space-time history of the string breakup, and with the presence of resonant hadronic states. 
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I. INTRODUCTION 

The Lund string fragmentation model [1] represents one 
of the most successful phenomenological descriptions of the 
non-perturbative stage of the hadronisation process, where 
partons (quarks and gluons) are converted into hadrons. The 
understanding of hadronisation is essential for a reliable in- 
terpretation of the hadronic data obtained in collider exper- 
iments, and the PYTHIA generator code J2l which contains 
the Monte-Carlo implementation of the Lund fragmentation 
model plays an important part in practically every analysis. 

The model uses the concept of a string with uniform energy 
density to model the confining colour field between partons 
carrying complementary colour charge. The string is viewed 
as being composed of straight pieces stretched between indi- 
vidual partons according to the colour flow. The fragmenta- 
tion of the string proceeds via the tunneling effect (creation 
of a quark-antiquarkpair from the vacuum) with a probability 
given by the fragmentation function. The space-time sequence 
of string break-up points defines the final set of hadrons, each 
built from a qq pair ( or a quark-diquark pair in the case of 
baryons) and a piece of string between the two adjacent string 
break-up points. The longitudinal hadron momenta stem di- 
rectly from the space-time difference between these vertices, 
see Fig.Q] 

The idea of a helix-like shaped string was first suggested in 
the study of the properties of soft gluon emission by Anders- 
son et al. JH. Under the assumption that the generating cur- 
rent has a tendency to emit as many soft gluons as possible, 
and due to the constraint imposed on the emission angle by 
helicity conservation, it was shown that the optimal packing 
of emitted gluons in the phase space corresponds to a helix- 
like ordered gluon chain. Such a structure of the colour field 
cannot be expressed through gluonic excitations of the string 
and it needs to be implemented as an intrinsic string property. 

The current paper is organized as follows: section II de- 
scribes the changes in the fragmentation model related to the 
helix string structure, the modification of the original helix 
string proposal, and the Monte-Carlo implementation of the 
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FIG. 1. Evolution of the QCD string in the rest frame of the qq 
pair. The endpoint partons loose their momentum as they separate 
and the string - the confining field - is created. The space-time 
coordinates of break-up vertices can be obtained from the relation 
[t,x]=(k + p++k~p~)/ic (k ~ 1 GeV/fm), with jfc+( _ J designating the 
fraction of parton momenta used for the build-up of the string field. 
The x direction is parallel to the string axis and the mass of partons 
is neglected. 



model. Section III uses a simplified version of the hadronisa- 
tion model for a detailed study of the experimental signature 
of the helix string model. Section IV deals with the charge 
asymmetry in the correlated hadron production. Section V 
provides an overview and discussion of the observed features. 
Appendices A and B show the results of first tuning studies 
using the LEP and LHC data. 



II. MODELLING OF TRANSVERSE MOMENTUM 

The implementation of a string with a helix structure rad- 
ically changes the way hadrons acquire their transverse mo- 
mentum. In the conventional Lund model (H], the transverse 
momentum of the hadron is the (vectorial) sum of the trans- 
verse momenta of the (di)quarks which were created via tun- 
neling during the breakup of the string. The transverse mo- 
menta of newly created partons are randomly sampled from 
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a gaussian distribution (with adjustable width) and their az- 
imuthal direction is uniformly random (i.e. two random num- 
bers are thrown at every string break-up point in order to gen- 
erate the associated transverse momentum, Fig. |2). 




FIG. 2. The generation of the transverse momentum of a hadron in 
the conventional Lund string model, obtained as a vectorial sum of 
the transverse momenta of quarks(diquarks) created in the tunneling 
process. The string breakup points are uncorrelated in the azimuthal 
angle and the (locally conserved) transverse momentum is sampled 
from a gaussian distribution. 

In the case of the helix ordered string, hadrons obtain their 
transverse momentum from the shape of the colour field itself, 
so that there is in principle no need to assign a momentum to 
new quarks in the string breakup. If we assume the gluon field 
has a form of an ideal helix and the string tension is tangen- 
tial to the helix, the transverse momentum stored in the string 
piece defined by two adjacent string break-up points can be 
written as 

p T = R sxp^ 71 ^ d<P (1) 
J<s>, 

where is the 'phase' of the helix (azimuthal angle) at 
the break-up point and R stands for the radius of the helix. 
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FIG. 3. In the helix string model, the transverse momentum is stored 
in the transverse structure of the gluon field. The tunneling process 
creates quark-antiquark pairs at rest and the string breakup points are 
correlated. 

The transverse momentum that the hadron carries is en- 
tirely defined by the properties of the helix field (described 
by two parameters: the radius and the winding density, see 
Fig. 0. The loss of the azimuthal degree of freedom in the 
string break-up is arguably the most significant consequence 
of the inclusion of the helix string model in the fragmentation 
process. 



ELI. Parametrisations of the helix string structure 

In the original helix string proposal [3], the phase difference 
of the helix winding is related to the rapidity difference of the 
emitting current by the formula: 

Ay 

A*=-p (2) 

where A<I> is the difference in helix phase between two 
points along the string, T is a parameter, and Ay is the rapidity 
difference which can be calculated as 

Ay = 0.51n(^), (3) 

with k + '~ representing the fractions of endpoint quark mo- 
menta defining a position along the string, see Fig[T] In this 
parametrization, the helix phase difference is related to the 
angular difference of points in the string diagram (FigJl^). 

An alternative parametrisation [4] sets the difference in the 
helix phase to be proportional to the energy stored in between 
two points along the string 

A<£ = (Ak + +AkT)M /2, (4) 

where Mq stands for the invariant mass of the string, 
^[rad/GeV] is the density of the helix winding, and Ak + = 
\kj — kj +l \, Ak~ — \kj — ^J+i I define the size of the string 
piece. This parametrisation describes a static helix structure 
where the helix phase remains constant at a given point along 
the string (FigHp). As discussed in Section 3.2 of [4], such a 
parametrisation emerges as a solution of equations developed 
in lU assuming the separation of soft gluons is dominated by 
the azimuthal component (neglected in J3t]). 

For practical reasons, only the static helix scenario (Eq. |4]) 
is studied in this paper - the original helix string proposal 
cannot be implemented for an arbitrary parton configuration 
without modifications ( the model contains a singularity at the 
gluon kink, where one of the fractions k becomes in Eq.[2]). 
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FIG. 4. The space-time evolution of the helix string structure in the 
rest frame of a qq string. Left: original proposal (Eq|2}, where helix 
phase evolves along the string axis in time. Right: modified proposal 
(EqlU, corresponds to a static helix structure. 
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For simplicity, the static helix scenario will be refered to as 
"helix string model" in the following, even though it repre- 
sents just one possible solution for the helix string parametri- 
sation. The phenomenology of the helix string depends on the 
exact definition of its shape and the experimental signatures 
may vary to some extent. 

II.2. Model implementation 

The helix string fragmentation model has been imple- 
mented in a PYTHIA-compatible form, and the fragmentation 
code of PYTHIA [50 has been adapted |6] to allow switching 
between the standard Lund string fragmentation and the he- 
lix string fragmentation (refered to as "standard" and "helix" 
fragmentation, respectively). 

The comparison of the helix string model with data requires 
the model to be extended to cover not only the simple case 
of a qq system but also an arbitrary multiparton configura- 
tion corresponding to the emission of hard gluons from the 
quark-antiquark dipole, as in Fig. [5] This is actually the most 
complicated part of the model implementation which requires 
some additional assumptions to be made. 

First, it is assumed that the helix phase runs smoothly over 
the gluon kink, i.e. the helix phases at the connecting ends of 
adjacent string pieces coincide. 

Second, it is assumed that the probability of the soft gluon 
emission depends on the mass of the string piece formed by 
colour partners in the quark-gluon cascade (the gluon momen- 
tum is equally split between the two adjacent string pieces), so 
that the helix phase difference between string endpoints, for a 
colour ordered system of N partons, becomes 

i<N 

A4> = J^£M i7 (5) 

where the sum runs over all (ordered) string pieces and M,- is 
the mass of the i-th string piece. In the simulation, the random 
choice of the helix phase at one point along the string defines 
the helix phase for the entire space-time history of the string 
evolution. In particular, the knowledge of the combined phase 
difference allows the fragmentation to proceed in the usual 
way, by steps taken randomly from one of the string's ends. 




FIG. 5. An illustration of the helix phase evolution of the modified 
helix model in the presence of a hard gluon kink on the string. 



III. EXPERIMENTAL SIGNATURE 

For a better understanding of the origin of observable ef- 
fects, it is useful to consider a simplified model configuration 
with ideal helix shape which provides a better illustration of 
the basic features of the model. The smearing of the helix 
shape can be obtained from the comparison with the data (Ap- 
pendix A). 

The correlations between transverse and longitudinal com- 
ponents of direct hadron momentum which originate from the 
helix string structure are best visible in a toy study of fragmen- 
tation of a simple qq string, shown in Fig. [6] The longitudinal 
direction and the transverse plane are defined with respect to 
the string axis. 




energy [GeV] 

FIG. 6. The dependence of the direct hadrons transverse momentum 
on their energy in the rest frame of the fragmenting helix string with 
radius R=036 GeV and ,y=0.7 rad/GeV. Toy model of hadronic Z° 
decay with suppressed parton showering. 



The helix string structure also creates a very distinct cor- 
relation between the azimuthal opening angle of two direct 




dE [GeV] 

FIG. 7. The dependence of the azimuthal opening angle of two 
direct hadrons on their energy-distance dE along the string in the 
rest frame of the fragmenting helix string with parameter y=0.7 
rad/GeV. Toy model of hadronic Z° decay with suppressed parton 
showering. For combination of hadrons with energy distance ex- 
ceeding a half-loop of the helix winding, two possible configurations 
exist. 
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hadrons and their distance along the string, which can be ex- 
pressed as a function of the amount of energy stored in the 
ordered hadron chain separating the pair, see Fig. [7] The po- 
sition of the hadron is associated with the middle point of the 
string piece which forms the hadron. 

Both types of correlations are rather difficult to observe ex- 
perimentally, though. In the presence of gluon kinks, it is im- 
possible to define the longitudinal axis of the whole string in a 
way which allows the separation of the intrinsic longitudinal 
and transverse momenta of all hadrons. On top of this im- 
portant source of smearing, it is necessary to account for the 
presence of non-direct hadrons (resonance decay products), or 
even - as in the case of hadron-hadron collisions - several over- 
lapping hadron systems, which contribute considerably to the 
dilution of the features associated with the underlying helix 
string structure. 

At this point, it is instructive to separate the effect of the 
parton showering (and string overlap ) from the effect of res- 
onance decay. The next step therefore consists in the study 
of a realistic configuration of input partons (quarks and hard 
gluons), but with fragmentation limited to the production of 
direct charged pions n~). 

Figure [8] shows the disappearance of the clear interference 
pattern between the transverse momentum and the energy of 
a direct hadron in a sample of hadronic Z° decays into light 
quark pairs (uu, dd), with parton shower enabled and reso- 
nance production disabled. The creation of heavier quarks in 
the parton shower and in the tunneling is suppressed, too. The 
difference between Figs. [6] and [8] is entirely attributed to the 
presence of a hard gluon cascade and the associated complex 
string topology. 
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FIG. 8. The dependence of the size of the transverse momentum(*) 
of the direct pions on their energy in the Z° rest frame. Hadronic Z° 
decay with suppressed resonance production. The underlying cor- 
relations are smeared by the parton shower.(*) transverse plane is 
defined with respect to the event Thrust axis. 

Even though the correlation pattern is smeared out, the 
presence of underlying correlations is still visible - indirectly - 
in the shape of the inclusive pj spectrum. Figure [9] shows the 
comparison of normalized p^ spectra for the standard frag- 
mentation and the helix string fragmentation, in the same sam- 
ple of Z°s decaying, ultimately, into a set of direct pions. The 
radius of the helix string has been adjusted (R = 0.5 GeV) in 



0.06 


1 1 


1 

- 


0.04 


Z° -> uD,d9 -> rm / 


\ 




standard fragmentation / 






helix fragmentation / 




0.02 






Pi 


i i 





1 0" 2 1 0" 1 1 

transverse momentum p [GeV] 



FIG. 9. The comparison of the inclusive transverse momentum (*) 
distribution obtained via standard fragmentation (full blue line) and 
via helix string fragmentation (dashed red line). Hadronic Z° decay 
with suppressed resonance production^*) transverse plane is defined 
with respect to the event Thrust axis. 



order to reproduce the average < pj > obtained in the con- 
ventional fragmentation (with a PT = 0.36 GeV). The parame- 
ter 5? is set to 0.7 rad/GeV and kept unchanged in this paper. 
The difference in the shape of the two distributions reflects the 
difference in model assumptions - the standard Lund fragmen- 
tation creates a gaussian-like shape while in the helix string 
fragmentation, the shape of the inclusive pj is closer to an 
exponentially falling function. The study of the shape of the 
inclusive pj spectrum is one possible way how to experimen- 
tally distinguish between the two model variants. Note that 
the difference between models is artificially enhanced in the 
sample where all hadrons are direct. 

A similar pattern is found in the study of 2-particle corre- 
lations, where the helix string structure causes an enhanced 
production of pairs of close direct hadrons. This is seen in 
Fig. [TO] in the distribution of the momentum difference Q de- 
fined as 




Q [GeV] 



FIG. 10. The comparison of the momentum difference distribu- 
tion for pairs of charged pions obtained via standard fragmentation 
(full blue line) and via helix string fragmentation (dashed red line). 
Hadronic Z° decay with suppressed resonance production. 
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Q = yj-(pi-pi) 2 , (6) 

where p,- stands for the 4-momenta of hadrons. The compar- 
ison is done using pairs of charged hadrons ( direct pions in 
our example ). 

The traces of the 2-particle correlation pattern shown in 
Fig.|7]can also be detected with the help of a suitably defined 
power spectrum on a selection of events with restricted parton 
shower activity. Indeed, this is what has been found in a recent 
study by ATLAS [7], where the sensitivity to the parameters 
of the helix string structure has been investigated on MC in 
low p7 inelastic proton-proton scattering. 

IV. CHARGE ASYMMETRY 

The prediction of low Q enhancement by the helix string 
model becomes even more intriguing when looking at the dif- 
ference between pairs of hadrons with like-sign and unlike- 
sign charge combination - the enhancement is significantly 
more pronounced for the unlike-sign pairs (Fig. [TTV 

O 1 1 1 ,— 



"O 




Q [GeV] 

FIG. 11. The comparison of the momentum difference distribution 
for pairs of charged pions obtained via standard fragmentation (blue 
lines) and via helix string fragmentation (red lines). The contribu- 
tions from pairs with identical charge ("like-sign") are indicated by 
dashed lines, the unlike-sign charge combinations correspond to full 
lines. Hadronic Z° decay with suppressed resonance production. 

A charge-asymmetric enhancement of 2-particle produc- 
tion at low Q is readily observed in hadronic data. However, 
this enhancement is seen in like-sign pairs and usually asso- 
ciated with the Bose-Einstein interference, though the phe- 
nomenon has never been successfully simulated from first 
principles, and a number of measurements actually contradict 
this interpretation (absence of correlations between overlap- 
ping hadronic systems H], JUL fToh . fiUl. d ependence of the 
size of the "source" on the particle type IU2I1 ). 

The presence of string-structure induced 2-particle correla- 
tions at low Q therefore raises questions : can these two ob- 
servations be connected ? is there a way to modify the charge 



flow in the model so that it reproduces the charge asymmetry 
in the data ? 

It turns out there are actually two ways how to influence the 
charge asymmetry of 2-particle correlations in the fragmenta- 
tion of the helix string, and both represent a natural extension 
of the helix string model, as explained below. 



IV. 1. Fluctuations in the space-time sequence of string 
breakup 

The PYTHIA generator fragments a QCD string from both 
of its endpoints, creating one direct hadron at each step, until 
the mass of the remaining piece of string falls below a cer- 
tain (adjustable) limit, and is then split into two final hadrons. 
This is a convenient way of implementing the fragmentation 
algorithm that neglects fluctuations in the lifetime of a string 
piece comprising several adjacent direct hadrons, as shown in 

Fig.na 




FIG. 12. The string area diagram showing the PYTHIA fragmenta- 
tion of qq string in an outside-in direction, creating a single hadron 
at each step (left diagram). The algorithm neglects fluctuations in 
the space-time sequence of the string breakup (right diagram) be- 
cause for the standard Lund string model, the final states resulting 
from both diagrams are indistinguishable in the energy-momentum 
space. The horizontal axis shows the distance along the string axis, 
the vertical axis indicates the time in the rest frame of the string. 

It is argued in [1], section 2, that neglecting these fluctu- 
ations does not affect the outcome of the modelling since 
the hadron final states are identical in the energy-momentum 
space. This argumentat is indeed valid for the standard 
Lund string model, where the transverse momentum of direct 
hadrons is acquired in the tunneling process and carried by 
constituent quarks (or diquarks). However, if we consider the 
helix string model, where the transverse momentum is car- 
ried by the string itself, the fluctuations do have an impact 
on the final hadronic system - the hadrons become effectively 
reordered along the helix gluon chain. This is illustrated in 
Fig. E] 

The impact of these fluctuations ( i.e. the presence of the 
string breakup configurations which are omitted in the stan- 
dard PYTHIA fragmentation chain ) on the relative contribu- 
tion from like-sign and unlike-sign pairs of hadrons is shown 
in Fig. [14] The ratio of like-sign and unlike-sign Q distribu- 
tions is plotted. Two scenarios are studied: the "enhanced" 
one, where the ordering is done in a way which favours the 
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FIG. 13. The string area diagram showing the effective swap of 
adjacent hadrons in the fragmentation of a QCD string with helical 
structure ( the vertical lines in the diagram show the evolution of a 
fixed helix phase). Since the transverse momentum of hadrons is as- 
sociated with the transverse structure of the string, the fluctuations of 
the space-time sequence of string breakup are reflected in the relative 
ordering of final hadrons in the azimuthal angle, creating distinguish- 
able hadronic final states. 



creation of "adjacent" like-sign pairs of hadrons, and the "nat- 
ural" one, where the fluctuations occur randomly according to 
the string area law |Q]|. In both scenarios, only the fluctua- 
tions relevant to the swap of adjacent (next-in-rank) hadrons 
are considered, and for the toy model of Z0 decay into direct 
pions, the average probability of the swap is around 0.4. For 
comparison, the ratios obtained with the standard PYTHIA 
ordering are included in the plot both for the standard Lund 
string and the helix string. 

IV.2. Resonances 

The presence of resonant hadronic states in the fragmen- 
tation chain has a significant impact on the size of correla- 
tions stemming from the underlying string structure and on 
the charge asymmetry. The principal question resides in the 
treatment of the resonance decay - the resonance, in particular 
a short-lived one, can be viewed as a piece of the string which 
preserves its properties, and decays according to string break- 
up rules, as a smooth continuation of the string fragmentation 
process. It is also possible that resonances, especially long- 
lived ones, gradually "lose" the memory of the QCD field. 

To see the impact of the different models of the reso- 
nance decay on 2-particle correlations associated with the he- 
lix structure of the string, another version of the toy modelling 
of hadronic Z° decay is set up, where the QCD strings are 
forced to decay into a set of direct p° resonances. The p°s 
consequently decay according to the usual PYTHIA decay al- 
gorithm ("std. p° decay"), or their decay obeys the helix string 



FIG. 14. The relative fraction of like-sign and unlike-sign pairs 
of charged pions as a function of their momentum difference Q, 
of the string type and of the ordering of hadrons along the string. 
The helix-string fragmentation favours the creation of close adja- 
cent unlike-sign hadron pairs in the conventional PYTHIA ordering 
(dashed lines). The asymmetry is attenuated in the presence of fluctu- 
ations which correspond to an effective reordering of hadrons along 
the gluon chain (full lines). Hadronic Z° decay with suppressed res- 
onance production. 



break-up rules ("helix p decay"). The comparison of the mo- 
mentum difference spectra for pairs of charged pions obtained 
with standard PYTHIA fragmentation and resonance decay, 
and with helix string fragmentation scenario combined with 
standard, resp. helix resonance decay, is shown in Fig. [15] It is 
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FIG. 15. The comparison of the momentum difference distribution 
for pairs of charged pions obtained via standard fragmentation and 
standard resonance decay (full blue line), via helix string fragmen- 
tation and standard resonance decay (dashed red line), and the he- 
lix string fragmentation incorporating the resonance decay (dotted- 
dashed green line). Hadronic Z° decay with forced p° resonance 
production. 



not too surprising to see that the hadronisation scenario where 
the resonance decays do not take into account the underlying 
string structure does not exhibit the sort of enhancement ob- 
served with direct pions - the final hadrons "lost" the memory 
of the generating QCD field in the decay of the intermedi- 
ate resonance state. It is also clear that the presence of res- 
onances has a dampening effect on the low Q enhancement: 
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FIG. 16. The ratio of the momentum difference of like-sign and 
unlike-sign pairs of charged pions obtained via standard fragmenta- 
tion and standard resonance decay (full blue line), via helix string 
fragmentation and standard resonance decay (dashed red line), and 
the helix string fragmentation incorporating the resonance decay 
(dotted-dashed green line). Hadronic Z° decay with forced p° reso- 
nance production. 



the offspring of a resonance decay has a momentum differ- 
ence Q 



Ml,. 



p Q — 4m^ so that the unlike-sign Q distribution 

contains a resonance structure, and the presence of a heavy 
resonance increases the average distance -along the string - be- 
tween other direct hadrons (compare Fig. [I5]and Fig. [ToV 

It is however quite interesting to see that the inclusion 
of resonances enhances the relative fraction of like-sign and 
unlike-sign pion pairs at low Q, as shown in Fig. [16] The 
effect is substantial and provides answers to questions formu- 
lated in the current section - the helix string model predicts a 
low Q enhancement dominated by like-sign pairs in the pres- 
ence of resonances which preserve the "memory" of the gluon 
field through their decay. 
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FIG. 17. The comparison of the momentum difference distribution 
for pairs of charged pions obtained via standard fragmentation and 
standard resonance decay (blue lines) and the helix string fragmen- 
tation incorporating the resonance decay (green lines). The full lines 
indicate contributions from unlike-sign pairs, the dotted lines corre- 
spond to the contributions from like-sign pairs of pions. Hadronic Z° 
decay with forced p° resonance production. 



Figure \T7\ compares the Q spectra, to stress the fact that 
even if the like-sign pair production dominates the low Q re- 
gion, some enhancement is also expected for the unlike-sign 
pairs when comparing the helix string scenario with the stan- 
dard Lund fragmentation. This feature may become instru- 
mental in the test of the model against the data. 



IV.3. Elongated source? 

Having successfully demonstrated - with the help of sim- 
plified fragmentation scenarios - that the enhanced produc- 
tion of like-sign hadron pairs with a low momentum differ- 
ence is compatible with the helix string model, the temptation 
is great to take a step further and to try to see if the model 
predicts some of non-trivial properties of these correlations. 
In the analysis of the LEP data fl3ll . it has been shown that 
the correlation function exhibits an asymmetry in the longitu- 
dinal and transverse direction - the phenomenon is known as 
'elongation' of the source in the terminology of Bose-Einstein 
interferometry. The asymmetry is seen as faster dampening of 
the correlation function along the longitudinal axis (the Thrust 
axis of the event). 




(GeV; 



FIG. 18. The comparison of the 2-dimensional momentum differ- 
ence of like-sign (top) and unlike-sign (bottom) pion pairs obtained 
with the helix string fragmentation model. Qi ons (Qside) correspond 
to the longitudinal(transverse) component of the momentum differ- 
ence. The circular band in the bottom spectrum is associated with 
the resonance decay. A clear enhancement of the Q s id e arm is visible 
in the like-sign spectrum. Hadronic Z° decay with disabled parton 
shower and forced p° resonance production. 
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Once more, the toy model of hadronic Z° decay with forced 
p° production and subsequent helix-like decay of p° reso- 
nance into charged pions (the model which exhibits significant 
enhancement of like-sign pair production at low Q) is used, 
but the parton shower is switched off to prevent the smearing 
of transverse and longitudinal components. The difference be- 
tween charged pion momenta is calculated in the longitudinal 
centre-of-mass frame of the particle pair (see [13] for defini- 
tion). Figure[T8]shows the 2-dimensional distributions of like- 
sign and unlike-sign combinations. The circular band which is 
visible in the unlike-sign spectrum corresponds to the p res- 
onance decay. The like-sign spectrum exhibits an enhanced 
particle production in low Qi ong region which is not visible 
in the unlike-sign spectrum, and which may be at the origin 
of the elongation effect (the detailed discussion of this fea- 
ture is left to a separate publication). Indeed, after subtraction 
of the two distributions, the resulting picture shows a certain 
'elongation' of the particle production, in any case an asym- 
metry in comparison of the shape of the like-sign pair excess 
in transverse and longitudinal directions (Fig. 119). 




[GeV] 



FIG. 19. The difference of the 2-dimensional Q distributions of 
like-sign and unlike-sign pion pairs obtained with the helix string 
fragmentation model exhibits an asymmetry which may be at origin 
of so-called "elongation" effect. Hadronic Z° decay with disabled 
parton shower and forced p° resonance production. 



V. OVERVIEW 

There are only a couple of free parameters in the helix string 
model: the helix radius R (which replaces the a Pr parameter 
in PYTHIA) and the density of the helix winding (parame- 
ter 5? in the case of the static helix scenario). The predictive 
power of the model has so far been used essentially for the 
development of the model itself : the necessity to accommo- 
date the charged asymmetry in the predicted 2-particle cor- 
relations leads to the reconsideration of "hidden" degrees of 
freedom of the string fragmentation model (the dependence 
on the space-time history of the string breakup). The exis- 
tence of these degrees of freedom translates into uncertainty 



of the model prediction which hopefully will be reduced by 
comparison with the experimental data. 

The key ingredient of the successful implementation of the 
helix model seems to be the adjustment of the relative fraction 
of direct hadrons, in particular short-lived resonances with de- 
cay conforming to the helix string fragmentation rules, and 
their ordering along the string. The suitable algorithm has yet 
to be developed, but some estimates of the observable effects 
can be made using a mixture of "ordinary" hadron production, 
and enhanced resonance production. 

An example of such an estimate is given in Fig.|20l for the 
study of the 2-particle correlations at LHC. The prediction of 
the helix string model is obtained from the combination of two 
samples: A (helix-string fragmentation without modification 
of the particle type content nor the resonance decay) and B ( 
helix-string fragmentation into p resonances followed by the 
resonance decay according to the helix string structure), in 
the proportion 4: 1 (the combined sample contains 20% of the 
enhanced resonance production). The proportion is chosen in 
such a way that the size of the predicted correlations does not 
exceed values observed in the LHC data 1 141. 



1.5 
1.4 
1.3 



PYTHIA8 minimum bias, non-diffractive 

— • — like-sign pairs 
— s — unlike-sign pairs 




Q [GeV] 

FIG. 20. An estimate of the contribution of the helix-string fragmen- 
tation model to the 2-particle correlation function in inelastic low p^ 
pp collisions at LHC, using PYTHIA8 generator with the helix-string 
option @]. 



For completeness, it should be mentioned that a signifi- 
cant improvement in the overall agreement of the MC predic- 
tions with the inclusive single particle distributions and event 
shapes measured at LEP has been obtained with the helix- 
string model in an earlier tuning study [!l5|. The study has 
been repeated with a Pythia8-compatible re-implementation 
of the model [6]. The results are reported in the Appendix A. 



V.l. Additional experimental input 

The density of the helix string winding can in principle be 
measured directly in the study of azimuthal correlations of 
charged hadrons. The experimental data are available H and 
the comparison with the predictions of the helix string model 
extended to the decay of short-lived resonances is shown in 
Appendix B. 
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Another experimental input, not discussed so far, may come 
from the study of the polarization of resonance decays. It fol- 
lows from previous sections that at least some of the short- 
lived resonances should preserve the helix string structure in- 
ternally and decay accordingly, and it should be possible to 
compare the helix-string model predictions with the measure- 
ments of the angular distributions in resonance decays. 

V.2. Additional theoretical input 

The concept of a helix-shaped confining gluon field of- 
fers a fresh insight into certain features observed in the data 
but not sufficiently understood from the theoretical point of 
view. The competition between the helix-string model and 
the Bose-Einstein interpretation of particle correlations is the 
most striking example, but there are other aspects of the prob- 
lem which merit a theorist's attention. The possibility of the 
existence of a helix-like ordered gluon field emerged from the 
study of optimal packing of soft gluons at the end of parton 
cascade, yet it might be that some ordered gluon pattern is al- 
ready present at the level of gluon emission from the leading 
quarks. If this is the case, the properties of the string, trans- 
mitted through the correlations of hadrons emerging from the 
interaction, offer a glimpse at the way the charged particle in- 
teracts with the surrounding field. 

On a more practical level, the study of the helix-string 
model would profit from further clarification of which helix 
pattern is preferred, or, if the theory does not have a prefer- 
ence, from the implementation of viable options in MC gen- 
erators so that the model variants can be confronted with the 
data. In case of the original helix-string model proposal [3], 
this implies the regularization of the singularity associated 
with the endpoint of the string and the extension of the model 
on multiparton string configurations. 

VI. CONCLUSIONS 

The phenomenological consequences of the replacement of 
the conventional pj modelling in the Lund string fragmenta- 
tion with an alternative model based on the helix structure of 
the QCD string are numerous and experimentally verifiable. 
It is shown how the constraints imposed by the helix-shaped 
gluon field translate into modifications of the shape of the in- 
clusive pj distribution and of the momentum difference of 
hadron pairs. The particle correlations predicted by the helix- 
string model can explain a large part of the enhanced produc- 
tion of close pairs of like-sign hadrons. The strength of the 
model consists in its high predictive power, due to a low num- 
ber of free parameters and the reduction of the number of de- 
grees of freedom in the modelling. 

APPENDIX A 

As a cross-check of an earlier tuning study lU5ll . the pre- 
dictions of Pythia 8-compatible implementation of the helix 



string model [6] are compared with the data measured by the 
DELPHI Collaboration B16I1 . The adjustment of model pa- 
rameters is done with help of Rivet lfl7ll and Professor IU8II 
packages. 

The study uses as input the inclusive particle spectra and 
the event shape variables as listed in II 1 511 . The model setup 
is identical to Pythia8 default (Tune:ee=3) except for fitted 
parameters shown in Table I, for the standard fragmentation 
scenario and the helix string fragmentation. 

The comparison of the fit results in the table confirms the 
conclusions reported in fll5ll . namely that the overall descrip- 
tion of the data improves when the helix string scenario is em- 
ployed, in particular for the subset of inclusive particle spec- 
tra. The fit provides an estimate of the smearing of the ideal 
helix shape (parameter HSF:sigmaHelixRadius, HSF stands 
for the HelixStringFragmentation class). The fitted values of 
parameters of the helical shape vary with the choice of the in- 
put data: the radius (HSF:helixRadius) between 0.4 and 0.54 
GeV, and the density of the helix winding (HSF:screwiness) 
between 0.54 and 0.91 rad/GeV. Similar spread of values have 
been obtained in 1 15] when studying difference between mod- 
elling using pj -ordered parton shower, and the ARIADNE 
parton shower. 



Tuned parameter 


Input data set 


Pythia 8 Tuned 


Event shapes & 


Inclusive particle 


(std.fragm.) 


incl.particle spectra 


spectra only 


StringPT:sigma 


0.276(1) 


0.264(3) 


StringZ:aLund 


0.315(4) 


0.262(9) 


StringZ:bLund 


0.689(6) 


0.60(2) 


TimeShower: alphaS value 


0.1418(1) 


0.1452(2) 


TimeShower:pTmin 


0.662(8) 


0.73(3) 


X 2 /N d0 f 


3.72 


3.87 


Pythia 8 


Event shapes & 


Inclusive particle 


+HELIX [6] 


incl.particle spectra 


spectra only 


HSF:screwiness 


0.918(4) 


0.54(1) 


HSF:helixRadius 


0.405(2) 


0.53(2) 


HSF: sigmaHelixRadius 


0.063(1) 


0.07(1) 


StringZ:aLund 


0.513(3) 


0.51(6) 


StringZ:bLund 


0.443(5) 


0.22(2) 


TimeShower: alphaS value 


0.1386(1) 


0.1382(4) 


TimeShower:pTmin 


0.767(5) 


0.74(3) 


X 2 /N d0 f 


2.93 


2.07 



TABLE I. Results of the tuning study using the DELPHI data fiUl . 



Among the inclusive particle spectra, the most significant 
improvement is obtained in the description of the mean trans- 
verse momentum measured as the function of the scaled mo- 
mentum (Fig|2l"li. the observable which is best suited to detect 
the correlations between the transverse and absolute momen- 
tum of hadrons. 
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X, 

FIG. 21. The comparison of the DELPHI data with tuned predic- 
tions of Pythia8 using the standard fragmentation (dashed blue line) 
and the helix string fragmentation (full red line). 



resonances, which means a significant increase of helix-string 
induced correlations in the model prediction. It is interesting 
to see how well these predictions compare with the measure- 
ment. 

The study is done in the following way: the Pythia 8 based 
implementation of the helix string model is used with the de- 
fault hadron mixture. The p and K* resonances decay ac- 
cording to their 'internal' helix string structure. The parame- 
ters describing parton shower,fragmentation function and the 
smearing of the helix radius are fixed to values obtained in 
the fit of Z° data (Appendix A). A selection of parameters de- 
scribing MPI, ISR and CR (see Table II) is retuned in order 
to achieve a good agreement between the data and the model 
prediction in the inclusive event selection and in the low-p^ 
depleted region, less sensitive to hadronisation effects. The 
main parameters of the helix structure are included in the fit 
and their optimized values are within the range obtained from 
the comparison of the model with Z° data. Since the data 
indicate a very small amount of colour reconnection, the re- 
connection probability is set to to stabilize the minimisation 
procedure. 



APPENDIX B 

Recently, the ATLAS Collaboration published the measure- 
ment of the azimuthal ordering of hadrons inspired by the 
search of the signature of the helix-like shaped QCD field [0] . 
The spectral analysis of correlations between the azimuthal 
opening angle and the longitudinal separation of hadrons pro- 
vides a possibility to measure the parameter describing the 
helix winding density in the model. The data show presence 
of correlations in the region of interest (known from the pre- 
vious study of hadronic Z° data). The interpretation of the 
observed effect is complicated by the fact that the discrep- 
ancies between the data and the predictions of conventional 
models cannot be attributed to the modelling of hadronisa- 
tion alone but require re-adjustment in the jet sector (sensitive 
to the amount of multiple parton interactions/MPI, ISR, and 
colour reconnection/CR). 

According to the arguments developed in this paper, the he- 
lix string model should be extended to the decay of short-lived 



Tuned parameter 


Input data set 


Pythia8 4C 


ATLAS.2012.I1091481 


+HELIX [6J 


inclusive + low-pr depleted 


HSF:screwiness 


0.61(2) 


HSF:helixRadius 


0.460(5) 


MultipartonInteractions:expPow 


3.7(1) 


Multipartonlnteractions :pT0Ref 


2.09(3) 


SpaceShower:pTORef 


1.81(4) 


BeamRemnants:reconnectRange 


0.[F] 


X 2 /N dof 


581/443 = 1.3 



TABLE II. Results of the tuning study using the ATLAS data |7] 
collected at t/s=1 TeV. 
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FIG. 22. The comparison of the ATLAS data with tuned predic- 
tions of Pythia8 using the helix string model extended to the decay 
of short-lived resonances (full line). To illustrate the effect of the 
helix string structure, the predictions of the model with the standard 
string fragmentation and resonance decay are indicated (dashed line), 
for the same parameter setup. Low-pr enhanced event selection (re- 
stricted by the cut on the max(p7-) < 1 GeV). 
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The predictions of the retuned model for the region most 
sensitive to the hadronisation effects (low-p^ enhanced re- 
gion) are shown in Fig |22] A reasonable agreement between 
the data and the model is seen in the oservable. The model 
somewhat underestimates the size of correlations in the 
power spectrum, and the predicted distribution shows a wavy 
structure. This is an indication that some of spectral com- 
ponents are missing in the simulation and it provides a hint 
for further improvement of the modelling (the origin of the 



discrepancy in the S-q may be related to a poorly adjusted dis- 
tribution of the invariant mass of hadronic systems). 

For comparison, the predictions of Pythia 8 with the same 
parameter setup but using the standard fragmentation and de- 
cay algorithm are shown together with the predictions of the 
helix string model. The difference serves as an illustration 
of the effect the helix string model has on the size of the az- 
imuthal ordering signal. 
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